The continuum description with pseudo-state wave functions 
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Benchmark calculations are performed aiming to test the use of two different pseudo-state bases 
on the Multiple Scattering expansion of the total Transition amplitude scattering framework. Calcu- 
lated differential cross sections for p- 6 He inelastic scattering at 717 MeV/u show a good agreement 
between the observables calculated in the two bases. This result gives extra confidence on the 
pseudo-state representation of continuum states to describe inelastic/breakup scattering. 
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Inelastic scattering at intermediate energies can be a 
useful tool to study multipole excitations of Borromean 
nuclei (like n Li and 6 He). Due to their loosely bound 
nature, to properly understand and interpret such re- 
actions, it is crucial to take into account the few-body 
degrees of freedom. At high energies, the Multiple Scat- 
tering expansion of the total Transition amplitude (MST) 
is a convenient framework that has already been applied 
to analyze such reactions for elastic 0, || as well as for 
inelastic [H, 0| scattering. In the latter case, the method 
can take into account spin excitations that occur when 
scattering from a spin target such as a proton. In these 
calculations, it is formal and numerically advantageous 
to represent the continuum states in terms of a basis of 
square-integrable functions, also known as pseudo-states 
(PS). Unlike the true scattering states, the PS vanish 
at large distances and hence the method will be only 
useful if the calculated observables are not sensitive to 
the asymptotic region. Moreover, calculations performed 
with different families of states, should converge to the 
same results, provide that enough states are included, 
and that the basis is complete within the radial region 
which is relevant for the process under study. 

Guided by this motivation, in this Brief Report we 
present benchmark calculations of proton inelastic scat- 
tering from 6 He within the MST scattering framework 
making use of two different PS bases to describe the 6 He 
continuum. We aim to check to what extent the calcu- 
lated breakup observables depend on the choice of the 
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PS functions. 

For a Borromean system, like 6 He, the wave function 
for a total angular momentum J (with projection M) 
and energy e, cp^ , can be expressed in terms of the 
Jacobi coordinates f (the relative coordinate between the 
valence nucleons) and R (the relative coordinate from the 
center of mass of the neutron pair to the core). 

It is also convenient to introduce a set of hyperspher- 
ical coordinates: the hyperradius p and five hyperspher- 
ical polar angles Q$ — {u,Q x ,<j) x ,9y,(f)y\. The former is 

defined as p = y/x 2 + y 2 with scaled coordinates x — 
2 _1 / 2 r and y — (2/V3)R. The angle a — arctan(ic/y) 
is the hyperangle and 9 x ,4> x ,0 y , <f> y the angles associated 
with the unit spatial vectors x and y. 

Within the PS method, the eigenstates ip^ M (r,R) are 
obtained by diagonalization of the Hamiltonian in a ba- 
sis of normalizable states. These states are conveniently 
expanded in a basis of Hyperspherical Harmonics of the 
form 

^(f,R) = R nl3 (p)Tj M (n 5 ) , (1) 
where Ti M (SI5) is the generalized angle-spin basis [j| 

r J M (n 5 ) = {y Ktx(yL (n 5 ) ® [ Xs2 ® x S3 ] s } JM (2) 

with Xsi the neutron spin functions and y^ci^i l(^s) the 
hyperspherical harmonics, 

yice x i y LM L (n 5 ) = ijiKUtv (a) [Yt, (x) ® Y ly (y)] ^ . (3) 

The functions il>Ki x l v have an explicit form in terms 
of Jacobi polynomials of the hyperangle a The set 
of quantum numbers — {JC£ x £ y LS} defines a channel, 
with £ x and £ y the orbital angular momenta associated 
with the Jacobi coordinates x and y, fC = £ x + £ y + 2u 
(i/=0,l,2, . . .) the hyperangular momentum, L = £ x + l y 
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the total orbital angular momentum and S the spin of the 
particles related by the coordinate x. In Eq. , R n p (p) 
are the hyperradial functions and n is an index that labels 
the basis states within a given channel (3. These functions 
are orthogonalized such that 



dp p 5 Rni3(p)R n > p{p) = <W 



(4) 



The aim of the present work is to compare two differ- 
ent choices for the functions R n p (p) in the calculation of 
breakup observables within the MST scattering frame- 
work. First, we consider the Gauss-Laguerre (GL) basis 
whose hyperradial part, R^ L {p), is given by 

R^ L {P) = Po- 3 [n!/(n + m l/2 Ll{z) exp(-z/2) , (5) 

with z — p/po, the generalized Laguerre polynomials, 
and po a parameter that sets the radial scale of the basis. 

The second choice is the Transformed Harmonic Os- 
cillator (THO) basis, recently introduced in Ref. [7[ for 
a three-body system. The THO method is based on the 
idea of transforming the bound ground state wave func- 
tion of the system into the ground state wave function 
of the Harmonic Oscillator (HO), defining a Local Scale 
Transformation (LST). The ground state wave function 
can be written as a linear combination of the basis func- 
tions |U), 



<p J ° Mo (r,R) 
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i?|(p)T^ M (Q 5 ) 
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where we have introduced the abbreviated notation 
tPo oM ° = ipi° Mo . Then, the equation that defines the 
LST for each channel /3 is 

~>')| 2 , (7) 
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where Rq® (s) is the hyperradial part of the HO ground 
state for the hyperangular momentum /C. Then, the 
THO basis is constructed for each channel applying the 
LST, s f3 (p), to the HO basis 
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(p)) 
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where 2 are generalized Laguerre polynomials of de- 
gree n. For channels not included in the ground state, 
information from one of the known (ground state) chan- 
nels with the closest quantum labels to the channel of 
interest is used to construct the LST, as explained in 
Ref. Q. 

Neither the GL nor the THO functions are eigenstates 
of the Hamiltonian, but they provide a complete and or- 
thonormal set in which the Hamiltonian can be diagonal- 
ized. For this purpose, the basis is truncated by setting 
a maximum value of the index n (n = 0, . . . , N) as well 
as a maximum hyperangular momentum /C max - Upon di- 
agonalization in the truncated basis, the eigenstates are 
obtained as 
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C^^(f,R). 
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where {ei} are their associated eigenvalues. 

From the derivation above, it becomes apparent that 
the GL basis is obtained in a more straightforward way 
than the THO basis. However, the latter has some ap- 
pealing properties that could make it more suitable in 
some situations. In particular, the THO basis has the 
advantage of being constructed from the ground state 
wave function of the system. Thus, when we diagonalize 
the Hamiltonian in a finite THO basis, the ground state 
is recovered for any size of the basis. By contrast, in the 
GL representation a large basis may be required to ob- 
tain a good description of the ground state. Also, note 
that the hyperradial part of the GL basis is the same 
for all the channels (3 while in the THO basis a different 
hyperradial part is calculated for each channel, with the 
correct behavior at the origin (p ). 

For a meaningful comparison between the two bases, 
we use the same three-body Hamiltonian to generate the 
GL and THO eigenstates for 6 He. In particular, we use 
the n-n potential of Gogny, Pires and Tourreil @ with 
spin-orbit and tensor components and we take the n- He 
potential from Q. Besides the pairwise interactions, an 
effective three-body potential is included, with matrix 
elements of the form 51 



U|» = 
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(10) 



l + (p/5)3 • 

The J=0 strength of this effective potential is tuned to 
reproduce the experimental three-body separation energy 
and the J > strength is adjusted to obtain the 
resonance at the experimental energy. 

We now consider the scattering process of 6 He, origi- 
nally in its ground state, |</?o° M °), *° a fi na l continuum 
state \ip^ M ), at excitation energy e and with total angu- 
lar momentum J (projection M), by means of its interac- 
tion with a proton, with initial (final) linear momentum 
k\ (k^) in the nucleon- nucleus center-of-mass frame and 
spin Si — 1/2 with projection a {&'). 

The double differential cross section for this process 
can be formally expressed as 
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where T denotes the transition amplitude operator [Iol |. 
This operator can be expressed as a multiple expansion 
series in the transition amplitudes tx for proton scatter- 
ing from each projectile sub-system X [ll|. At high ener- 
gies and for small momentum transfers, this expansion is 
expected to converge quickly. If only the leading term of 
the series is retained, the single scattering approximation 
(SSA) is obtained [1 03: 
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(12) 
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with 1=2,3 for the halo neutrons, and X=4 for the core. 
The proton - X subsystem transition amplitude satisfies 
the Lippmann-Schwinger equation 



(13) 



with v\x the interaction between the nucleon and I sub- 
system. Within the impulse approximation, the propa- 
gator Go = (E + — K) 1 contains the kinetic energy op- 
erators of the proton and all the projectile subsystems. 

Here E is the kinetic energy, E 
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in the overall 

center of mass frame, and fiNA is the proton-projectile 
reduced mass. 

Within the PS method, the scattering states ip J t M in 
Eq. are approximated by the pseudo-states <p^ M . 
Hence, the double differential cross section (fTTj) becomes 
a single differential cross section for each pseudo-state, 
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where we have replaced the T matrix operator by its sin- 
gle scattering approximation. Making use of the impulse 
approximation [2| , the matrix elements for the scattering 
for each constituent can be further simplified, leading to 
the following factorized form for the scattering from one 
valence nucleon (Z=2): 

E( TO3 t> ?Tl4 
t[b0S p s> p ][ui2, A) x p W -,s T su,\ I — A >— — A 



1,0 



2:; 



M; 
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(15) 



with M23 — m 2 + TO3, M234 = m 2 + m 3 + to 4 and where 
we have introduced the momentum transfer A = k ^ — ki 
and the energy parameter 0J12 0] and where S p = {Si<j} 
(S' p — {Si<j'}) are the incoming (outgoing) spin of the 
nucleon and its projection and, St = {Jo-Mo} (S T = 
{JM}) the initial (final) total spin of the halo valence 
pair. The amplitude t[bps p s^] is given in terms of the 
tensor components of the nucleon-nucleon transition am- 
plitude @, [l3| . The transition density form factors, 
P[b(3;S T s' T ei] j depend exclusively on the structure of the 
composite system. Its explicit expression as a function 
of the hyperradial parts of the wave functions of the ini- 
tial and final states, can be found in 0. 

The scattering from the core, assumed here as spinless, 
can equivalently be written as: 



t[00S p S' p ] (k>14, A) X P[00;S t St £ >] 
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(16) 
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FIG. 1: (Color online) Angular differential cross section for 
the breakup of 6 He on protons at 717 MeV per nucleon, lead- 
ing to J n = + continuum states of the 6 He nucleus. The 
three lines represent the calculation with the THO basis, for 
several values of N (indicated by the labels). Each curve in- 
cludes the contribution from eigeinstates up to 10 MeV in 
excitation energy, according to Eq. (|17[l. 



where, as before, W14 is the appropriate energy parameter 
0. The angular differential cross section for 6 He inelas- 
tic scattering (breakup) is then obtained by summing all 
excited states contributions, 



V- d(7 JJo 

^ dn 



(17) 



For the evaluation of Eqs. ([15]) and ([16]) one needs the 
(free) transition amplitudes for proton scattering from 
the valence nucleons and the core. For the former, we 
used the NN Paris interaction. The transition amplitude 
for the a core was generated from a phenomenological 
optical potential, of Woods Saxon form, with parameters 
obtained by fitting existing data for the elastic scatter- 
ing ofp+ 4 He at E p = 700 and 800 MeV, as detailed in 
Ref. [|. 

We first study the convergence of the breakup observ- 
ables with respect to the basis size. For this purpose, we 
consider the THO basis, truncated at different values of 
N. The maximum hyperangular momentum was set to 
A^max = 20. This yields the three-body force parameters 
y 3B = -2.4 MeV, for J = 0, and Vp = -0.85 MeV, for 
J>0. 

In Fig. [T] we show the angular distribution of the cal- 
culated inelastic differential cross sections. For simplic- 
ity, only the J" = + continuum is included, and the 
Coulomb interaction between the proton and the a core 
is ignored. The three lines represent the SSA calcula- 
tion for different values of the basis size, according to the 
choice of the parameter N. The three cases are in al- 
most perfect agreement, indicating that in this reaction 
the convergence with the basis size is very fast. 

Next, we study the dependence of the breakup observ- 
ables on the choice of the basis, by comparing the cal- 
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FIG. 2: (Color online) Calculated contributions for breakup 
differential cross section leading to final states with J 71 — 
+ , 1~, 1 + and 2 + , using the SSA approximation. Solid and 
dashed lines refer to the GL and THO basis, respectively. 

culations in the GL and THO representations. As be- 
fore, the maximum hyperangular momentum was set to 
K max = 20, and only eigenstates below 10 MeV are con- 
sidered. The index n was truncated to N — 20 and N = 4 
for the GL and THO bases, respectively. With this model 
space, the number of pseudo-states in the GL (THO) ba- 
sis is: 31 (30) for 0+, 63 (86) for 1~, 53 (49) for 1+ and 79 
(81) for 2 + . For the GL basis, the range parameter was 
set to pa — 0.25 fm, which provides a basis that extends 
up to about 20 fm in the hyperradius. With these pa- 
rameters, the ground state obtained after diagonalization 
of the Hamiltonian appears at -0.9781 MeV and -0.9549 
MeV, for the GL and THO bases, respectively. 

In Fig. [5] we compare the inelastic angular distribu- 
tions calculated in the GL and THO bases. The separate 
contributions for J^=0 + , l - , 1 + and 2 + final states are 
also shown. As before, the Coulomb interaction is ne- 
glected. The thick lines are the incoherent sum of all 
these J 71 ' contributions. Solid and dashed lines corre- 
spond, respectively, to the calculations with the GL and 



THO bases. For each curve, the contribution of eigen- 
states up to e max = 10 MeV are added incoherently, ac- 
cording to Eq. (|T7|) . We consider only the forward angles 
@c.m. < 30° since the SSA is not expected to work well 
at large momentum transfers [2]. We see that, at these 
angles, the dominant contribution to the breakup cross 
section comes from the l - states, while for 8 c . m , > 25°, 
the 2 + excitation becomes dominant. Finally, the popu- 
lation of the unnatural 1 + states is almost negligible at all 
angles. We notice that this excitation mode requires spin- 
flip transitions which, according to these calculations, are 
very small in this reaction. 

For all these contributions, the GL and THO bases pro- 
vide very similar results, suggesting that the calculated 
observables do not depend on the choice of the continuum 
representation, provided that enough states are included. 

In summary, in this Brief Report we have calculated 
proton inelastic scattering from 6 He at 717 MeV/u, us- 
ing as scattering framework the single-scattering approx- 
imation and two different pseudo-state representations of 
the 6 He continuum: the GL and the THO. Provided that 
enough states are included, both bases predict essentially 
the same inelastic differential cross section. Furthermore, 
the studied observables converge very quickly with the 
size of the basis. These results support the reliability 
of the pseudo-state method as a useful and convenient 
tool to treat scattering problems dealing with continuum 
states. This analysis could be extended to other PS bases 
and reactions. Furthermore, it could be applied to other 
scattering frameworks, for which the PS method has been 
also implemented, such as the continuum discretized cou- 
pled channels (CDCC) method [MODI. 
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